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Abstract 

Employing the standard hard-scattering approach and the running couphng method we 
calculate a class of power-suppressed corrections ~ l/Q^", n = 1,2,3... to the electromagnetic 
71*^7 transition form factor (FF) Q'^F^^^^Q'^) arising from the end-point a; — > 0, 1 integration 
regions. In the investigation we use a hard-scattering amplitude of the subprocess 7 -|- 7* q + Q, 
symmetrized under the exchange fij^ /I^ important for exclusive processes containing two 
external photons. In the computations the pion model distribution amplitudes (DA's) with 
one and two nonasymptotic terms are employed. The obtained predictions are compared with 
the CLEO data and constraints on the DA parameters 62 (/^o) ^^'^ ^4(^0) normalization 
point fiQ = 1 GeV^ are extracted. Further restrictions on the pion DA's are deduced from the 
experimental data on the electromagnetic FF Ft^{Q'^). 



PACS numbers: 12.38.Bx, 13.40.Gp, 14.40.Aq 



I. INTRODUCTION 



The tt" meson electromagnetic transition form factor (FF) F^^^IQ"^) is among the simplest 
exclusive processes for investigation of which at large momentum transfer the perturbative 
QCD (PQCD) methods Q i, 3 

can be applied. Because of the recent CLEO data j^, 
where the form factor F^^^Q"^) was measured with high precision, the interest in this process 
has been renewed. Thus during the last few years for computation of F^^(Q^) the various 
theoretical methods and schemes were proposed [1, H, 0, H 0|- The aim here is twofold: 
to elaborate methods for the calculation of F^^jlQ"^) within PQCD and, at the same time, 
to extract from experimental data information on the pion distribution amplitude (DA). 
The latter, being independent of a specific exclusive process and universal quantity, is an 
important input ingredient in studying various processes that involve the pion. 

It is known that in experiments the 7r°7 transition was explored at momentum trans- 
fers of ~ 1 — 10 GeV^, which are far from the asymptotic limit oo, where the 
PQCD factorization formula with the pion leading- twist asymptotic DA leads to reliable 
predictions. In the present experimentally accessible energy regimes, power-suppressed cor- 
rections ~ n = 1,2,.. play an important role in explaining the experimental data 

1^. There are numerous sources of power corrections to F.,^y{Q^). For example, the pion 
higher-twist (HT) DA's and higher Fock states generate such corrections. Power corrections 
can also originate from the end-point regions a; — *• 0, 1 as a result of the integration of the 
PQCD factorization expression with the QCD running coupling a^iQ^x) [as{Q'^x), x = 1—x] 
over the pion's quark longitudinal momentum fraction x. In fact, in order to reduce the 
higher-order corrections to a physical quantity and improve the convergence of the corre- 
sponding perturbation series, the renormalization scale //^ {V^'r), i-e., the argument of the 
QCD coupling, in a Feynman diagram should be set equal to the virtual parton's squared 
four- momentum ^Ifl]. In the exclusive processes the scale (/I/j) chosen this way inevitably 
depends on the longitudinal momentum fractions carried by the hadron constituents. For 
the photon-meson transition we have = Q^x and Jlj^ = Q^x, because at two leading order 
diagrams of the partonic subprocess 7 + 7* ^ q + q, absolute values of the virtual quark 
and antiquark squared four-momenta are determined by these expressions. But then the 
PQCD factorization formula diverges, since ag(Q^x) [as{Q^x)] suffers from end-point x ^ 
[x — > 1] singularity. The running coupling (RC) method solves this problem by using a 
Borel transformation and applying the principal value prescription. As a result, one obtains 
the Borel resummed expression for the 7r°7 transition FF, which contains power-suppressed 
corrections. The RC method in conjunction with the infrared (IR) renormalon calculus^was 
used for computation of such power corrections to the 7r°7 and 177, r]''-/ transition FF's 

to the electromagnetic FF's of the light mesons Fm (Qi) (M = 71, K, pl) mis, he 

well as to the gluon-gluon-r/' meson vertex function '15]. 

In the present work we compute power corrections to the 7r°7 transition FF employing the 
version of the hard-scattering amplitude symmetrized under replacement yu|j ^ 7l|.. The 
symmetrization procedure is important for exclusive processes with two external photons 
(gluons) in the hard-scattering Feynman diagrams, because it allows one to treat within 
the RC method both virtual and real photons (gluons) on the same footing. The latter is 
required in order to consider the 7r°7* and 7r°7 transitions in a unifying way, i.e., to get in 
the limits — > 0; 1 {u is the asymmetry parameter) from the n^'y* transition FF F^^. (Q^, u) 
the FF Fj^^lQ"^) of the 7r°7 transition. The advocated method was used in our previous work 
jisj l to investigate the virtual and on-shell gluon-77' meson transitions. In what follows we 
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refer to this approach as the symmetrized RC (SRC) method. 

This paper is organized as follows: in Sec. II we introduce the symmetrization procedure 
of the hard-scattering amplitude for the 7r°7 transition. Here we present our results for 
the Borel resummed [Q^F^^(Q^)]'''^* FF obtained within the SRC method. Section III is 
devoted to detailed analysis of the Q"^ oo limit of [Q"^ Ft^^^^Q"^)]^^^ . In Sec. IV we perform 
numerical computations and from comparison of our predictions with the CLEO data extract 
constraints on the parameters 62(1 GeV^) and 62(1 GeV^), 64(1 GeV^) in the pion DA's with 
one and two nonasymptotic terms, respectively. Further restrictions on the DA's arising from 
analysis of the pion electromagnetic FF are described in Sec. V. In Sec. VI we make our 
concluding remarks. 



II. THE PHOTON-MESON TRANSITION FORM FACTOR 

A. The symmetrized version of the hard-scattering amplitude 

The real photon-pseudoscalar M meson electromagnetic transition FF FM-y{Q'^) can be 
defined in terms of the amplitude F^"^, 

T'^-' = te^FM^{Q'')e^'""'Paqip, (2.1) 

for the process 

7*(gi)+ 7(^2) -^M(P), (2.2) 

where = —qf is the momentum transfer. 

At the large momentum transfer the FF -Fm7(Q^) is given by the factorization formula 
of the standard hard-scattering approach (HSA) 

Fm^{Q^) = [Th {x, Q\ + T| [x, fil)] ® 0m(x, /i|). (2.3) 

Here the function Th{x, Q^, /i|^), 

Th (x, Q', fil) = Th (x, Q', fil) + T| (x, Q^ fil) , (2.4) 

is the hard-scattering amphtude of the subprocess 7 + 7* — g + g , 0m(x, fj.p) is the meson 
DA, //|, is the factorization scale, and x = 1 — x, x being the longitudinal momentum fraction 
carrying by the meson's quark. In Eq. (j2.3|) the shorthand notation 

Th (x, (5^ (S) 0m(x, fip) = Th (x, Q^, fip'j 4>m{x, ^p)dx (2.5) 

is used. 

It is evident that a physical quantity, represented by the factorization formula, Eq. ()2.3|) 
being a sample one, does not depend on renormalization and factorization schemes and scales 
employed for its calculation. But at any finite order of the QCD perturbation theory, due 
to truncation of the corresponding perturbation series, the hard-scattering amplitude (j2.4p 
depends on both the factorization and renormalization /i^ scales. Since higher-order 
corrections in PQCD computations, as a rule, are large for both inclusive and exclusive 
processes, in order to get reliable theoretical predictions within the PQCD by means of 
the truncated perturbation series, an optimal choice for these scales, i.e., a choice that 
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minimizes higher-order corrections, is required. The factorization scale /zl^ in exclusive 
processes is traditionally set equal to the momentum transfer Q^, because higher-order 
corrections contain terms ~ \n{fi'jp / Q"^) and such choice eliminates them in hard-scattering 
amplitudes. As a result, in the factorization formula only a hadron DA explicitly depends 
on the scale fij^ = Q"^. 

The situation with the renormalization scale /x^ is more subtle. Really, this scale appears, 
in general, not only in higher-order corrections to the hard-scattering amplitude, but also 
determines the scale of the QCD coupling as(/i|j). In order to reduce higher-order corrections 
to a physical quantity, in exclusive proceses the scale should be taken equal to the square 
of the momentum transfer carried by a virtual parton in each leading order Feynman diagram 
of the underlying hard-scattering subprocess [lO|. For the real photon- meson transition these 
scales are determined by the leading order diagrams of the subprocess 7 + 7* — g -|- g and 
are given by the expressions 

fil = Q\ Jll = Q'x. (2.6) 

After these remarks let us turn to our formulas (j2.3p and (j2.4j) . In accordance with the 
"tradition," in this work we set = and in what follows omit the dependence of the 
hard-scattering amplitude on the scale u%- Then, for the hard-scattering amplitude at the 
next-to-leading order (NLO) we get [1^ 



N 1 
Q^x 



l + Ci 



a. 



Ml 



tix) 



where the function t{x) is given by the expression 



xlnx 



t{x) = In a; — 



9. 



X 



(2.7) 



(2.8) 



Here N is the constant, which depends on the quark structure of the meson, Cp = 4/3 is 
the color factor. The second function T^(a;, Q^, yu|.) can be obtained from Eq. ()2.7|1 using 
the replacement x ^x 

T2(x,g^/i^)=Ti,(x,Q^ 71^). (2.9) 
The hard-scattering amphtude Th{x,Q^, fij^) must be symmetric under exchange x 



(2.10) 



The replacement x ^ x,hj means of which the function Tfj{x, Q^, /i|.) is found, in general, 
has to be applied also to the renormalization scale changing it to Jl'j^. In the standard 
HSA one treats the /ifj and 7l|. scales on the same footing by setting them equal, as a rule, 
to . The choice = /I^ = satisfies both requirements ()2.9p and ()2.10p important 
for the hard-scattering amplitude. In the framework of the RC method the scales and 
71^ have to be chosen in accordance with Eq. ()2.(j|l . Then the function Tjj{x, Q^, /i^ = Q^x) 
takes the following form 



N 1 
Q^x 



1 + C^^tix) 
An 



The second part of the hard-scattering amplitude is given by the expression 



(2.11) 



TUx,Q' 



N 1 



l + C,^^tix) 
An 



(2.12) 
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One can see that within the RC method the requirements ()2.9|) and ()2.1U|) hold as welL 

In the framework of both the standard HSA and RC method the M7 transition FF can 
be calculated employing the formula 

= 2T^(x,Q2)®<^^(a;,Q2). (2.13) 

In the last step we take into account that the DA of the pion is a symmetric (pAii^jQ"^) = 
(pM^x^Q"^) function. 

The M = 7r°, 77, 77' meson electromagnetic transition FF's were computed within the 
RC method in Refs. 0|- In this work we generalize our approach by performing the 
computation of the transition FF in the context of the RC method, but employing 
instead of Eqs. (|2.11|) and (j2.12|) their versions symmetrized under /i^ ^ Jij^ exchange, i.e., 

Thix,Q') = {1 [«,(g2a;) +a,(g%)] t(x)| , (2.14) 

and 

Tl (x, = (i + k(g'x) + «s(g'x)l t(x)] . (2.15) 



Q^xi Sir 

In the standard HSA Eqs. ()2.14|) and ()2.15p coincide with Eq. ()2.7p and its x ^ x partner 
T^(x, g^,/i|. = g^), respectively. It is also not difficult to demonstrate that Eqs. fl2.9p . 
(I2.1()|l . and ()2.13j) hold for the hard-scattering amplitude determined by the new functions 
Tjj{x,Q^) and T2(x,g2) . 

Here some comments concerning the symmetrization procedure are in order. To clarify 
this important point let us note that the virtual and real photons enter into the considering 
process fl2.2j) in an unequal manner. Indeed, the M7 transition FF FM-yiQ"^) depends only 
on g^ = —qf (^2 = 0). At the same time the virtual photon- meson, 

7*(9i)+7*(g2)^Af(P), 

transition FF FM-y* {Q"^ , uj) is a function of the photon total virtuality and asymmetry 
parameter u (see the second paper of Ref. [3|), 



Q' = Ql + Ql c. = |' 



2 ■ 



In the limits uj 0;1 the equality 

FMr(g',c^ = 0;l) = FM-,(g') 

must be valid. In order to meet this requirement and describe the real and virtual photon- 
meson transitions within the RC method in a unifying way, we adopt in this work Eqs. ()2.14|l 
and ()2.15p . because in the limits a; — >• 0; 1 the form factor Fm^j* {Q"^ , uj) found in the context 
of the RC method leads to FM'yiQ^), computed by means namely of Eqs. (I2.14|) and ()2.15p . 
The symmetrization procedure, being discussed here, was used in Ref. ^] to calculate the 
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virtual and on-shell gluon-?]' meson vertex function. In the present work we concentrate on 
the FF FM-yiQ'^), leaving the detailed analysis of FM'y*{Q'^,uj) for a future publication. 



B. The pion distribution amplitude 



Calculation of the FF F^^(Q^) requires the knowledge of the pion DA ^^^(x, Q^), which 
is one of the key components in Eq. ()2.3|) . It is known ^3] that the pion DA can be 
expanded over the eigenfunctions of the one-loop Brodsky-Lepage equation, i.e., in terms of 
the Gegenbauer polynomials {C^/^(2x — 1)}, 



n=2A.. 



where (j)asy{x) is the pion asymptotic DA, 

4>asy{x) = V3fnX{l - 

with = 0.0923 GeV being the pion decay constant. 



X) 



(2.16) 



(2.17) 



The evolution of the DA on the factorization scale Q is governed by the functions bn{Q 



as(/io) 



00 



(2.18) 



In Eq. ()2.18|) {7n} are anomalous dimensions defined by the expression 



7n = Cf 



(n + l)(n + 2)+ fr!,j 



(2.19) 



The constants &n(/io) = b^ are input parameters that form the shape of DA's and can be 
extracted from experimental data or obtained from the nonperturbative QCD computa- 
tions at the normalization point /Iq. The QCD coupling constant as(Q^) at the two-loop 
approximation are given by the expression 



«s(g' 



An 



/3oln(QVA2) 



1 - 



2(3i lnln(QVA2 



(2.20) 



Here A is the QCD scale parameter, Pq and Pi are the QCD beta function one- and two-loop 
coefficients, respectively. 



Po 



11 



-n 



Pi 



51 Hf. 

3 ^ 



In the limit Q"^ ^ oo all model DA's (|2.16|) reduce to the asymptotic form (pasyi^)- The 
nonasymptotic terms ~ C^/^(2a; — 1), n > 2 determine the deviation of the pion DA from the 
asymptotic form at moderate energy regimes and depend on the nonperturbative mesonic 
binding effects. 

For the pion in the literature the various phenomenological DA's were proposed (gLItI. [TtI 
isl . 19l |. Thus, for example, in Ref. [l^, employing the QCD sum rules method, the following 



pion DA was predicted 



X, fil) = (pasy{x) 1 + 0.758C2^/^(2a; - 1) + 0.3942C4^^(2x - 1 



3/2/ 



(2.21) 
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where the normahzation poin is /io = 0.5 GeV. 

The coefficients 6° ^4 were also extracted from the CLEO data on the 7r°7 transition 
FF in Ref. 01 • The authors used the QCD hght-cone sum rules approach and included 
into their analysis the NLO perturbative and twist-four corrections. They found that in the 
model with two nonasymptotic terms, at the scale /io = 2.4 GeV, the pion DA has the form 

= (pasyix) [l + 0.19C2^^{2x - 1) - 0.14cf ^(2a; - 1)] . (2.22) 

As is seen the pion DA's extracted from the experimental data depend on the used 
methods and on their accuracy. Although one claims that the meson DA is a process- 
independent quantity describing the internal structure of the meson itself, exploration of 
different exclusive processes with the same meson leads to a variety of DA's. This means 
that employed methods have shortcomings or do not encompass all mechanisms important 
for a given process. Such situation is pronounced in the case of the pion. The investigation 
carrying out in this work intends to improve the situation with the n^'j transition FF by 
taking into account at least one class of power corrections to the FF 

To proceed it is convenient to expand the DA ()2.16p over x and rewrite it in the following 
form: 

oo 

Q') = 4>asyix) ir„x", (2.23) 
n=0 

where the sum runs over all n. The new coefficients Kn in the case of DA's with two 
nonasymptotic terms are given by the expressions 

i^o = 1 + 662(Q') + 1564(Q'), = -SOMQ^) + 7hiQ% 

K2 = 30[62(Q') + 28bi{Q% K3 = -60 ■ 2164(Q'), = 30 ■ 21b^{Q^). (2.24) 

Here the functions &2(<5^) and b4^{Q^) are defined by Eq. ()2.18p with 'J2/P0 and ■Ja/ Po being 
equal to 

72 ^ 50 74 ^ 364 

/3o 81' f3o 405' ' 



and 

72 2 74 364 



= 4, 



Po 3' Po 375' 
below and above the charm quark production threshold, respectively. 



C. The vr 7 transition FF within RC method 

Computation of the photon-pion transition FF F^^((5^) implies, naturally, integration 
over X in accordance with Eq. ()2.13|1 . Having inserted the explicit expression of the hard- 
scattering amplitude T^(x,(5^) ()2.14j) and the pion DA ()2.23j) into Eq. ()2.13j) we encounter 
divergences, arising from the singularities of the coupling constant as{Q'^x) and as{Q^x) in 
the limits x ^ 0; 1. In the standard HSA this problem is solved by freezing the argument of 
the coupling constant and performing corresponding integrations with as(Q^) [or as{Q'^/2)]. 
In the RC method we allow the QCD coupling to run and therefore have to propose some 
method to cure these divergences. 
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As the first step we express the running coupling as(Q^x)^, in terms of Os^fQ^). This aim 
can be achieved by applying the renormalization-group equation to asiQ'^x) [20]. As a result 
we find 



as(Q2)/?iln[l + lna;/t] 



(2.25) 



271(30 1 + lnx/t 

where a^iQ"^) is the one-loop QCD coupling constant and t = 47r//3oas(<5^) = ln(Q^/A^). 
Equation ()2.25|) expresses as{Q'^x) in terms of as(<5^) with an ~ a'^{Q'^) order accuracy. 

Inserting Eq. ()2.25p into the formula for the transition FF Eq. ()2.13p . we obtain integrals, 
which are still divergent, but can be calculated using existing methods. One of them (see 
for details Ref. allows one to obtain the form factor as a perturbative series in as(Q^) 
with factorially growing coefficients C„ ~ (n — 1)!, 



n=l 



Air 



(2.26) 



But, it is known that a perturbative QCD series with factorially growing coefficients is 
a signal for the IR renormalon nature of the divergences in Eq. (j2.26|) . The convergence 
radius of such series is zero and its resummation should be performed by employing the 
Borel integral technique. Namely, one has to determine the Borel transform B [Q'^Ft^^] (u) 



of the corresponding series j21| 



B 



7Vy 



[u] 



E 

n=l 



U 



n-l 



n 



1)! 



Crii 



(2.27) 



and in order to define the sum Eq. ()2.26p . or to find the resummed expression for the form 
factor, one has to invert B [Q^F^^] {u) to get 



Q Fj^'yi^Q 



P.V. 



duexp 



B 



Q'^F, 



TT'J 



(2.28) 



Because the coefficients of the series Eq. (I2.26|) behave like C„ ~ (n — 1)!, the Borel transform 
(|2.27p contains poles located at the positive u axis of the Borel plane, which are exactly the IR 
renormalon poles. Therefore the inverse Borel transformation (j2.28|) can be computed only 
after regularization of these pole singularities. One of the methods of such regularization, 
adopted also in the present work, is the principal value prescription. In other words, the 
IR renormalon divergences in Eq. ()2.28p have to be removed by computing the integral in 
the sense of the Cauchy principal value. Only after this regularization the inverse Borel 
transformation defines the resummed FF. 

Fortunately, these intermediate operations can be omitted with the help of the following 
operations. Namely, let us introduce the inverse Laplace transformations of the functions in 
Eq. ^rm . i.e.. 



1 



{t + zY 



1 f°° 

— — / duexp[—u{t + zYu"'^ , Reu > 0, 
z/ Jo 



(2.29) 



and 



^ Similar consideration is valid also for the runnig coupling as{Q^x) 
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\n[t + z] 

JtTW 



POO 

/ (iwexpf— u(t + z)](l — 7^ — lnM)u, 
Jo 



(2.30) 



where r(^) is the Gamma function, je — 0.577216 is the Euler constant, and z = Inx [or 
z = \Rxm the case of as{Q^x)]. Then, using Eqs. ()2.29j) and ()2.3()j] for the QCD coupling 
OisiQ^x) we find 



Ajr fOO 

a,{Q^x) = — / rfMe""*i?(M,t)x-", 

PO "'0 

where the function R{u, t) is defined as 

2 Pi 

R{u, t) = 1 ^"^^(1 — 7s — Int — In-u). 

Po 

Having used Eq. ()2.3H) and performed integration over x, for the scaled 7r°7 transition 
FF we get 



(2.31) 



(2.32) 



^ ( K 4 /-oo 



u,t) 



Y^Kn \An{u)+An{u) 



n=0 



(2.33) 



Here the term ~ An{u) appears in the result of the integration of the second term in Eq. 
()2.14|) . whereas the term ~ An{u) owing to the third term in Eq. ()2.14p . The functions 
An{u) and An{u) have the following forms: 



An{u) 



(n 



and 



AM 



u 



+ 



1 
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(n + 2-M)3 (n + 2-M)2 (n+l-u)(n + 2-u)' 
:B{2 - u, (3) - ^B{1 - u, (3) \p=n+2 - 95(2 - M, n + 1) 



(2.34) 



= B{n + 1, 2 - m) [(^(n + 1) - -^{n + 3 - uyf + ip'in + 1) - ip'{n + 3 - u) 

- B{n + 2,1- u) [ip{n + 2)-ip{n + 3- u)] - 95(2 -u,n + l), (2.35) 

where B{x, y) is the beta function -B(x, y) = T{x)T{y) /T{x + y) and ip{z) = d[\n T{z)]/dz. 

The functions An{u) and An{u) contain the poles on the positive real axis of the plane u. 
Indeed, the function A„(m) has the finite number of triple, double, and single poles located 
at the points Uq = n + 2 and triple and single ones at Uq = n + 1. In order to reveal the pole 
structure of the function An{u), it is convenient to use the following formulas [23^ 



oo -j^ 

^(.)^-7. + (.-l)E^^^^^. 



(k + zy 



(2.36) 
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Here we write down, as an example, the function Aq{u), which after some manipulations 
takes the form 



Ao(u) = (2-u) 



E 



.k= 



-y t 



2 



1^1 1 - 8 

Now it is clear that Aq{u) contains the infinite number of the double poles located at 
uq =^ k + 3 and the single ones at = 1, 2, /c + 3. The similar analysis can be fulfilled 
for An{u), n > as well. Hence by employing Eq. ()2.3H) we have transformed the end- 
point X —>■ 0; 1 divergences in Eq. ()2.13j) into the IR renormalon pole divergences in Eq. 
(I2.33|) . The integral in Eq. (|2.33|) is the inverse Borel transformation (|2.28p . where the Borel 
transform of the NLO part of the scaled FF is defined (up to constant factor) 

as 

oo 

B[Q^F^^] (u) ~ Riu, t) J2 Kn [An{u) + An{u)] . (2.38) 

n=0 

The IR renormalon divergences in Eq. ()2.33|) must be removed by means of the the principal 
value prescription. The inverse Borel transformation after such regularization, as we have 
just pointed out above, becomes the resummed form factor [Q'^Ft,^{Q^)Y^'^. Therefore all 
integrals over u hereafter have to be understood in the sense of the Cauchy principal value. 



The expression [Q'^ Ft^^{Q^)Y^^ contains the power-suppressed corrections ~ 1/Q 



2n 



n 



1,2, .. . to the scaled FF, implicitly existing in the QCD factorization formula (j2.3j) . The 
detailed discussion of relevant problems can be found in Refs. [l^ . Here, for completeness, 
we outline the important points of this analysis. To make the discussion of this question as 
transparent as possible, let us for a moment neglect the nonleading term ~ al{Q'^) in Eq. 
()2.25|) and consequently make the replacement R{u, t) — 1 in Eq. ()2.31|) . Then the integrals 
in the scaled and resummed FF with multiple IR renormalon poles at Mq = ^ can be easily 
expressed in terms of the integrals with a single IR renormalon pole at the same point [see 
Eqs. (|3.6|) and (j3.15|) ]. so that our formula (|2.33p will consist of some linear combinations of 
the integrals, 

' 72n(Q), (2.39) 



j3o Jo n — u n" 
where f2n{Q) are the moment integrals, 

UQ) = ^jykkP-'a,{e). (2.40) 
The integrals fp{Q) were calculated in Ref. ji^] using the IR matching scheme: 



fp{Q) 



N 



n=0 



2'Kp 



as(Q') 



[n!-r(n + l, pln(g//i,))], (2.41) 



where /i/ is the infrared matching scale and r(n+l, z) is the incomplete Gamma function. In 
Eq. (j2.41|) {/p(yU/)} are phenomenological parameters, which represent the weighted average 



10 



of as(fc^) over the IR region < k < fij and act at the same time as infrared regulators 
of the right-hand side (RHS) of Eq. ()2.39p . The first term on the right-hand side of Eq. 
fl2.41|) is the power-suppressed contribution to fp{Q) and models the "soft" part of the 
moment integral. It cannot be calculated within the perturbative QCD, whereas the second 
term on the RHS of Eq. (|2.41|) is the perturbatively calculable part of the function fp{Q), 
representing its hard perturbative "tail." In other words, the IR matching scheme allows 
one to estimate power corrections to the moment integrals by explicitly dissecting them 
out from the full expression, and introducing new nonperturbative parameters fp{fii). The 
same moment integrals fp{Q) computed in the framework of the RC method [LHS of Eq. 
()2.39|) ]. contain information on both their soft and the perturbative components. Therefore 
we can state that the scaled and resummed FF ()2.33|) contain power corrections ~ l/Q^"- 
In phenomenological applications both the IR matching scheme and the RC method can be 
employed. But the RC method has an advantage over the IR matching scheme, because 
it allows one to compute the functions fp{Q) without introducing the new nonperturbative 
parameters /x/ and fp{fii). Moreover, using this method, the parameters fp{fii) themselves 
can be calculated in good agreement with model calculations and available experimental 
data [iai24j]. 

But the principal value prescription itself generates in the each integral over u higher-twist 
ambiguities, 

where $q(Q^) is a calculable function fixed by the residue of the integral at the pole uq = q 
and Nq is some numerical constant. The ambiguities taken into account in Eq. ()2.33|) modify 
the Borel resummed 7r°7 transition FF, yielding 

[Q'F^,{Q')Y^^ ^ [Q'F^,{Q')Y^' + [Q'F^,{Q')f''. (2.42) 

The HT term depends on the known functions {^q{Q^)} and coefficients {Kg} and on 
the unknown numerical constants {Ng}. In accordance with the "ultraviolet dominance 
assumption" this HT ambiguity allows one to estimate higher-twist corrections to the scaled 
form factor Q"^ Ft^^{Q'^) coming from sources another than end-point integration. By fitting 
the constants {A^^} to experimental data one can deduce some information concerning the 
magnitude of such corrections. 



III. ASYMPTOTIC LIMIT OF THE RESUMMED 7r°7 TRANSITION FF 

As we have emphasized above, the resummed 7r°7 transition FF contains the power 
corrections appearing due to the end-point integration. These corrections in the region of 
a moderate momentum transfer are essential for explaining the experimental data |^ 
T^ . 13| . But it is also evident that in the asymptotic limit ^ oo , where all higher- twist 
corrections should vanish, the standard HSA with frozen as{Q'^) and the pion asymptotic 
DA (l)asy{x) leads to the correct expression for the 7r°7 transition FF. Consequently, in the 
limit Q"^ ^ oo from the resummed FF we have to regain the asymptotic one. 

In the limit Q"^ ^ oo the pion DA goes to its asymptotic form, i.e., 

Q^)*^-^ (l)asy{x), (3.1) 
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which in the terms of the coefficients Kn means 

Kq^I, Kn^O, n> 0. (3.2) 

We take also into account that in this hmit the subleading term in the expansion of as{Q'^x) 
through as(Q^) has to be neglected. In other words, in the limit — oo we have to fulffil 
the replacement 

/ due-'''R{u,t)^^ due-""'. (3.3) 
Jo Jo 

After these operations the resummed FF takes the following form: 

[Q^F^^{Q')Y'' "^"-^ v^/.iV |l + ^J^ due~^' [M^) + M^)] } • (3-4) 

But Eq. fl3.4p is not the final expression, because in the integral above t = ln((5^/A^) and 
its — > cxo limit still has to be computed. 

To this end, we start from the simple case and consider the integral 



/■oo /■ 

Ji = / due-^'Aoiu) = 
Jo Jo 



oo 

ut 



due' 



1 9 9 

+ 



(1-m)3 (2-m)3 (2-m)2 1-u 2-u 

(3.5) 

The integrals in Eq. ()3.5|) with the triple and double poles can be reduced to ones with 
single poles: 

roo e-'^^du 1 In A \n^Xh{X'^) 



(n-uY 2n2 2n 2 A" 
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10 {n — uY n A"- Jo n — u A" 

where the ffist two equalities are obtained performing integrations by parts. Here the loga- 
rithmic integral is defined as 

H(O^PyJ^;^^, (3.7) 

and A = Q^/A^. Now using the expansion of li{C"')/C"' inverse powers of In^ jl5l |. 

C" nlnC ^0 ('^I'^O 
and keeping in the expressions 

,U{X-) ,^MX-) 



In'A— InA- 



A" ' A" ' 

terms up to 0(1/ In A) order, we get 

The situation with the second integral, 

/■oo ^ 

/2 = / c/Me-"*Ao(M), (3.10) 
Jo 
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is more subtle. In this case, instead of using the exphcit form of Aq{u), we consider the 
integral 

due-""' / dxx^~''t{x) = / due-""' / dxx^-%x), (3.11) 
JO Jo Jo 

from which Eq. (jH.lOj) has been derived. We are going to explain our technique, analyzing 
one of the components of the function t{x). Namely, let us calculate the Q"^ ^ oo limit of 
the integral 

Il= due-""' rfa;x^~"ln^(l -x). (3.12) 
Jo Jo 

Having expanded ln^(l — x) in powers of x, 

oo -| 

-x) = 2Y: T—rmk + 1) - ^(l)]x^+\ (3.13) 
fc=i + 

we obtain 

nl OO 1 1 

k=l 

Substituting Eq. (j3.14j) into the integral Jj, 

^^2g— (3.15) 
and using the leading order term in expansion (|3.8|) . in the limit — > oo we get 



»1 oo 1 1 

I fai-'ln=(l-x)=2i:— |#.(fc+l)-<;'(l)l^-j:33^. (3.14) 



(3.16) 



Now, having repeated the described above operations in the reverse order, it is easy to see 
that 

oo 1 1 1 1 /•! 



In other words, the asymptotic limit Q"^ ^ oo transforms the integral in accordance with 
the rule 

jiQ^ /'c^a;xln2(l-x). (3.18) 



In A JO 

The same conclusion is valid also for the other terms from Eq. ()2.8|) . Summing up, we derive 
the limit of the integral I2, 

h ^ I' dxxt{l -x) = j^i?(l, 2) { [(^(1) - ^(3))2 + ^'(1) - ^'(3); 

-[^(2) -^(3)] -9}. (3.19) 

The expression ()3.19p without the factor 1/lnA is nothing more than Aq{u) at m = 0. The 
following operations are trivial and lead to 

13 



As is seen, both the functions /i and I2 have the same hmits. Consequently, for the resummed 
FF we find 



in deriving of which the value of the constant A^, 

N = 



1 - ^asiQ') 

on 



12 et - e 



(3.21) 



(3.22) 



for the pion has been utilized. Equation ()3.21|) can be readily obtained within the standard 
HSA employing the pion asymptotic DA. Our analysis proves that in the asymptotic limit 
the Borel resummed FF leads to the correct expression ()3.2ip . which we consider as one of 
justifications of the symmetrization procedure. It is also worth remarking that the "old" 
version of the hard-scattering amplitude Th-(x, Q^) [see Eqs. fl2.11|l and ()2.12j) ] gives the 
correct asymptotic FF as well, that is evident from Eq. ()3.9|1 . Hence in the asymptotic 
limit both the ordinary and symmetrized RC methods describe correctly the n^'y transition 
FF, the difference between them being sizeable at the moderate values of the momentum 
transfer, ~a few GeV^. 



IV. EXTRACTING THE PION DA FROM THE CLEO DATA 

In this section we present the pion phenomenological DA's extracted from the CLEO 
data within the SRC method. In our calculations below we shall use the following values of 
the parameters A and /io 

A4 = 0.25 GeV, = 1 GeV^ (4.1) 

As is known (see for review Ref. jl^), the IR renormalon calculus can be applied for 
the estimation of power corrections to some physical quantity in the region of the high 
momentum transfers ^ A^. Our choice for the parameters ()4.H) leads to the requirement 
16 ^ 1. Because the recent CLEO data 01 on the 7r°7 transition FF lie in the domain 
1.64 GeV^ < < 10 GeV^, we include them into our numerical analysis to deduce the pion 
model DA's. Namely at these moderate momentum transfers the power corrections play the 
important role, modifying both qualitatively and qualitatively predictions for Q'^F.^.y^Q^) 
obtained within the standard HSA. 

The Borel resummed n^'j transition FF implies summations over n and k, the latter 
arising from ~ An{u) terms. The summation over n does not create problems, because in 
our studies we use the pion asymptotic and model DA's with one and two nonasymptotic 
terms. Therefore the maximal value of n in the sum is Amax = and 2, 4, respectively. It is 
worth noting that Eq. ()2.33p is a general expression valid for the pion DA's with an arbitrary 
number of nonasymptotic terms. The next terms ~ C„(2a; — 1), n > 4 can be easily included 
into our scheme by modifying only expressions of the coefficients Kn and A^max- Contrary 
to the case with n, at fixed n summation over k runs from k = to k = 00 and has to 
be truncated at some /cmax- In other words, the results of numerical computations depend 
on fcmax- In order to check their sensitivity to a chosen value of fcmax? we have performed 
calculation of the FF with fcmax = 50 and fcmax = 100. We have found that for the pion 
asymptotic DA the ratio 

o^^2^ ^ [g^F.,(g^)]'-^-(fcmax = 100) 

^ [Q'F^,{Q^)Y-^{k^^ = 50) ^ - ' 
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FIG. 1: The scaled and resummed transition form factor Q"^ Ft^^{Q'^) as a function of . 
The curves ASY and ASY* are computed using the pion asymptotic DA ()2.17|) . The curve ASY* 
corresponds to the FF obtained within the ordinary RC approach, whereas for calculation of the 
curve ASY the SRC method is employed. The upper dashed_ line shows the model-independent 
Q"^ ^ oo limit for the FF. The data are borrowed from Ref. 



at = iGeV^ is equal to GeV^) = 1.0027 and to /^(lOGeV^) = 1.0011 at = 
10 GeV^. We have obtained a similar picture employing the pion various model DA's. Since 
the correction to Q'^F^^^IQ'^) originating from the next k = 51 — 100 terms does not exceed 
3 ■ 10~^ of those from the first k = — 50 ones, in numerical computations we set fcmax = 50. 
Such output is understandable, because in the resummed FF dominate contributions from 
the nearest to m = IR renormalon poles. 

We start our analysis of the 7r°7 transition FF from the pion asymptotic DA in order to 
reveal the impact of the symmetrization procedure on the predictions, as well as to find out 
how large is the deviation of these predictions from the data points. Our results are shown in 
Fig.lU As is seen, the scaled and resummed FF Q^-FV^(Q^) computed using the SRC method 
in the region of the momentum transfers 1.64 GeV^ < < 5 GeV^ are considerably larger 
than the one obtained by means of the ordinary RC approach. As a result, the deviation 
of the curve ASY from the data points are smaller than that of ASY*. Nevertheless, such 
deviation exists and some admixture of nonasymptotic terms in the pion DA is needed to 
explain the data. 

In Fig.|21we depict (the shaded area) the la area for values of the input parameters 6°) i>1 
in the 6° — 63 plane. This means that the n^'y transition FF computed in the context of the 
SRC method employing the pion model DA's with Gegenbauer coefficients belonging to the 
shaded region describes the CLEO data with a la accuracy. 

In Fig. El we plot the la area for the 7r°7 transition FF itself. The central curves with 
62 = 0.16, 64 = (the DA's with one nonasymptotic term) and 62 = 0.25, 64 = —0.05 (the 
DA's with two nonasymptotic terms) are also shown. One sees that the shapes of the curves 
with 64 = and 64 < differ from each other. Indeed, the curves with 6^ < are sharper 
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FIG. 2: The la areas in the 64 — 62 plane of the input parameters. The shaded area is found from 
analysis of the CLEO data on the ir^^ transition FF. The region bounded by the dashed Hnes is 
extracted from the data on the pion electromagnetic (EM) FF. The dot-dashed line is the diagonal 
determined by Eq. (|4.4|) . 
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FIG. 3: The scaled and resummed ir^'y transition FF vs Q^. The shaded area demonstrates la 
region for the FF. Correspondence between the curves and the input parameters is; for the solid 
line ^2 = 0.25, 64 = —0.05; for the dashed line 62 = 0.16, 64 = 0; and for the dot-dashed line 
= 0.23, b1 = -0.05. 
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FIG. 4: The scaled and resummed form factor Q'^Ft^^{Q'^) with and without HT ambiguities. The 
sohd hue describes FF without HT ambiguities. The broken hues are found employing Eq. ()2.42|1 
and numerical constants {Nq} shown in the figure. 

relative to ones with 6° = 0. Therefore the boundaries of the la area are determined by 
superposition of curves of these two types. Our analysis in the case of the DA's with one 
nonasymptotic term leads to the following estimation 

6^ = 0.16 ±0.02, 6^ = 0. (4.3) 

In the case of DA's with two nonasymptotic terms allowed valus of fog? ^4 cover the shaded 
area in Fig. |21 Below we write down sample values of the parameters, 

6° = 0.2, 6° G [-0.02, -0.04], 
bl = 0.25, 6° G [-0.05, -0.07], 

and 

bl = 0.3, b^ G [-0.09, -0.1]. 
The "diagonal" of the la area is determined by the expression 

bl + 1.466^ = 0.16, bl G [0, -0.11] (4.4) 

and is shown in Fig. El by the dot-dashed line. 

As we have noted above, the principal value prescription generates HT ambiguities that 
in conjunction with the "ultraviolet dominance assumption" can be used to estimate HT 
corrections to the form factor originating from another source (for example, from the pion 
HT DA's). We have performed relevant computations, as a sample, for the pion DA with one 
nonasymptotic term 6°(1 GeV^) = 0.16 and {Ng} = ±0.6, g = 1, 2, . . . , 50 (Fig.i)). We find 
that the values {Ng} = ±0.6 determine the upper and lower bounds for the constants {Ng} 
in order that FF remain within the la region. The 7r°7 transition FF with HT ambiguities 
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corresponding to {Ng} = —0.6 {{Ng} = 0.6) at < 3.5 GeV^ is larger (smaller) than the FF 
without such corrections and are smaller (larger) for > 3.5 GeV^. The HT ambiguities, 
obeying the "la constraint" do not exceed the level ~ ±5% of the transition FF at the 
momentum transfers = 1.64 — 2 GeV^ and reach only ~ =pl.8% at = 9 — 10 GeV^. 



V. THE ELECTROMAGNETIC FORM FACTOR F^{Q'^) 

In this section we compute the pion electromagnetic FF F-,^{Cp) in the framework of the 
RC method in order to extract further constraints on the parameters h^, h\. 

The FF Ft^[Q'^) is the important quantity characterizing the pion, which was thoroughly 
investigated in the context of PQCD jl^, [iaH^]. It was also studied in various experiments 



28l |29|. Within the RC method FF's Fm(Q^) of the light mesons M = K, pi were 
considered in Refs. 0, 0, 0,1131 • Therefore below we outline only main stages of the RC 
analysis of the FF Fm(Q^). 

In the standard HSA the FF F.„^{Q'^) is given by the factorization formula 

F^iQ') = dx dycPliy, Q^)Th {x, y, Q') 0.(x, Q^). (5.1) 

Here Th{x, y, Q"^) is the hard-scattering amplitude of the subprocess gg'+j* — > qq', = —q^ 
is the momentum transfer, q being the four-momentum of the virtual photon. In Eq. (|5.ip 
the factorization scale from the very beginning is chosen equal to p\ = Q'^. 
At the leading order of PQCD the amplitude Th{x, y, Q^) has the form 



TH{x,y,Q ) = 



2 a(/^|) ^ 1 a(^|; 



xy 3 xy 



(5.2) 



In accordance with the ideology of the RC method the argument of the QCD coupling in 
Eq. ()5.2p has to be chosen as 

p% = xyQ^, JI% = xyQ^. (5.3) 

Such choice allows one to get rid of terms ~ \n{xyQ^ /jlj^), In^xyQ"^ / fij^) appearing in the 
amplitude Th{x, y, Q"^) at the next-to- leading order of PQCD and minimizes the higher-order 
corrections to Ft^[Q'^). We can also adopt the scheme 

p\ = xQ^, JIr = xQ^, (5.4) 

obtained from Eq. ()5.3p by freezing y. In the framework of the standard HSA one freezes 
both of X, y and compute the form factor with pji = ]I% = (or Q^/4 ) [26^ . 

In the above we have chosen x as the running variable. Alternatively, we can fix x and 
choose y as the running one or to compute the mean value of the sum of the FF's calculated 
using both possibilities; due to the symmetry of Th{x, y, Q"^) and Eq. 1)5.11) itself with respect 
to X, y we will get the same result. Of course, the second option ()5.4|1 leaves in the NLO 
correction some logarithmic terms, but it leads to better agreement with the experimental 
data than the choice ()5.3p [13^ . Therefore in our computations we use the option ()5.4|) . Since 
in the considering process the partonic hard subprocess contains only one external virtual 
photon, we do not perform the symmetrization of the hard-scattering amplitude. Stated 
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differently, below we write down the expression [Q"^ Ft^{Q'^)Y'^'^ , obtained in the framework of 
the ordinary RC method. 

The Borel resummed pion electromagnetic FF are determined by the formula [Ic 

[Q'F^m]"'' = £ KiB{2 + /, 1) 



POO 

y^Kn due-''^R{u,t)B{2 + n,l-u). (5.5) 



oo 

X 

n=0 



The integrand in Eq. ()5.5p has a finite number of the single IR renormalon poles. In fact, 
the sums in the general expression ()5.5|) in practice run up to some Lmax, A^max, which for 
DA's with two nonasymptotic terms are Lmax = A'^max = 4. The maximum number of IR 
renormalon poles results from the term ~ -8(6, 1 — u). The latter can be rewritten in the 
following way 

r(6)r(i-^) ^ 120 

^' ^ r(7-^i) (l-n)(2-M)...(6-M)' 

making our statement evident. It is implied that the pole divergences are removed by the 
principal value prescription. 

In the asymptotic limit —>■ oo from the resummed FF we recover the standard HSA 
expression. In fact, acting along the line described in the detailed form in Sec. Ill, we get 

[Q2F^(Q2)]re. Q!^°o (1M^5(2, 1) r ciMe-"*5(2, 1-u). (5.6) 

From Eqs. ()5.6|) and ()3.8p we obtain 

IQ^F^mr^^'^^ IGnf^a.m, (5.7) 

which can be found in the context of the standard HSA by employing the pion asymptotic 
DA. 

To perform the numerical analysis of the scaled and resummed pion FF Q'^F^^IQ'^) and 
extract constraints on the pion DA's from such consideration, we need to specify the exper- 
imental data that will be used in the fitting procedure. Unlike the 7r°7 transition FF, where 
we have precise CLEO data for large momentum transfers, the situation with the Q^F^^iCp) 
is somewhat controversial. Thus the corresponding data were obtained indirectly from the 
pion electroproduction experiments through a model-dependent extrapolation to the pion 
pole. Moreover, the points > 2 GeV^ are imprecise suffering from the large errors and, in 
addition, there are big gaps between data points themselves. The data on the FF Q^-F^(Q^) 
reported recently by the Collaboration do not change the whole picture, because the 
highest value of at which the measurements were performed is = 1.6 GeV^. There- 
fore, to improve the precision of the la analysis under the circumstances, we include into 
our fitting procedure data points > 1.18 GeV^ and slightly exceed in this way the range 
of validity of the RC method. But because our curves describe the data at such low values 
of as well, we find our approach justified. Finally, let us note that the datum point 

= 9.77 GeV^ is also included into our scheme and it strongly restricts the la region. 

The results of our numerical calculations are plotted in Figs. and |21 The la area for 
the pion scaled electromagnetic FF and the central curve with the Gegenbauer coefficients 
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Q'(GeVl 

FIG. 5: The pion scaled and resummed electromagnetic FF Q'^Ft^{Q'^) as a function of . The 
shaded area is the lu region for the form factor. The data are taken from Refs. 28j| (the circles) 
and [2^ (the rectangles). In the la analysis only the solid data points are used. For the central 
solid line the input parameters are 62 = 0.23, 64 = —0.05. For comparison the FF obtained by 
means of the asymptotic DA is also plotted. 

62 = 0.23, 64 = —0.05 are demonstrated in Fig. El The la region for the parameters h^, h\ 
of the pion DA's is shown in Fig. |21 They obey, for example, the following constraints 

hi = 0.16, bl e [-0.045, -0.05], 

6° = 0.2, 6° G [-0.039, -0.058], 

and 

bl = 0.28, bl G [-0.047, -0.061]. 

The overlap of the la regions in Fig. |21 determines the la area in the plane 64 — 62, within 
which both the ir^j transition and the pion electromagnetic FF's are in agreement with the 
corresponding data at the level of a la accuracy. As is seen, this area is rather restricted 
and the values of the parameters 6° ^4 

6^ = 0.235 ±0.035, 6° = -0.05 ^ 0.01. (5.8) 

The FF Q'^Fj^lQ'^) is more sensitive to HT ambiguities than the tt^'j transition FF. Actu- 
ally, in Fig. inifor 62 = 0.23 and 64 = —0.05 the scaled FF, corrected by the HT ambiguities, 
is plotted. These ambiguities for {Ng} = 0.1 and {Ng} = —0.08 reach ±(5.6 — 3.6)% of the 
FF in the region ^ 12 - 1.6 GeV^ and ±1% in the domain 9 - 10 GeV^ It is worth 
noting that the estimation of the HT ambiguities are obtained within the "la constraint" 
(see the previous section). In general, admissible values of the constants {Ng} and of the 
input parameters b^, 64 are strongly correlated. 
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FIG. 6: The pion electromagnetic FF with HT ambiguities. The sohd hne is the original FF with 
62 = 0.23, 64 = —0.05. The broken lines include the HT ambiguities with constants {Nq} = 0.1 
(the dashed line) and {Nq} = —0.08, (the dot-dashed line) g = 1, 2, ... 6. 

VI. CONCLUDING REMARKS 

In this work we have calculated the power corrections to the 7r°7 transition FF, originat- 
ing from the end-point regions x — > 0, 1 due to integration of the standard HSA factorization 
formula with the QCD running coupling over the longitudinal momentum fraction x, car- 
rying by the pion's quark. To this end, we have employed the RC method combined with 
techniques of the IR renormalon calculus. We have used the symmetrized under replacement 
A^R 71^ version of the hard-scattering amplitude of the partonic subprocess 7* + 7 — > 9 + ^- 

We have obtained the Borel resummed expression [Q"^ F^^{Q'^)Y'^'^ for the transition FF. 
For this purpose in the inverse Borel transformation we have removed IR renormalon di- 
vergences by means of the principal value prescription. Each IR renormalon pole uq = n 
in the Borel transform B^^"^ Ft^^]{u) corresponds to power correction ~ l/Q^" contained in 
the scaled and resummed FF. Since, in the considering process the Borel transform has an 
infinite number of IR renormalon poles, the expression ()2.33p . in general, contains power 
corrections ~ l/Q^"", n = 1,2, ...00. In numerical computations we have truncated the cor- 
responding series at nmax = 50. As an important consistency check, we have proved that the 
result obtained within the SRC method in the asymptotic limit — > 00 reproduce the stan- 
dard HSA prediction for the transition FF. This provides justification for the symmetrization 
procedure applied in the RC method. 

We have compared our predictions with the CLEO data and obtained restrictions on 
the input parameters of the pion DA's with one and two nonasymptotic terms. Further 
constraints on the admissible set of DA's have been extracted from the data on the pion 
electromagnetic FF F-j^lQ"^). We have concluded that the pion DA's with the parameters 
Eq. ()5.8|) describe the experimental data on both the Q'^ Ft,j{Q'^) and Q'^Ft,{Q'^) FF's with 
the la accuracy. 



21 



It is important that DA's extracted from the CLEO data are suitable for explanation 
of the pion electromagnetic FF, whereas in the context of the standard HSA for describing 
of these FF's one has to pose on the pion DA contradictory restrictions (or to model soft 
contributions to using mechanisms beyond the scope of the perturbative QCD). In 

fact, in the framework of the standard HSA the pion asymptotic DA considerably under- 
estimates the data on the electromagnetic FF In order to cover a gap between 
the data and theoretical curves one has to introduce model DA's with positive and large 
input parameters, the Chernyak-Zhitnitsky DA being one of the prominent examples. 
At the same time, (pasy{x) overestimates the CLEO data on Q'^Ft,^{Q^) and, on the contrary, 
model DA's with negative input parameters are needed. The RC method solves this problem 
due to power corrections taken into account in both of these quantities. Really, the power 
corrections arising from the end-point integration regions at moderate momentum transfers 
significantly enhance the pion electromagnetic FF 0, They also enhance the absolute 



value of the NLO contribution to the FF F^^^IQ'^). Since the contribution of the NLO term 
to F^^ylQ"^) is negative, power corrections effectively reduce the leading order contribution 
to FF. It turns out that for some model DA's these effects lead to a satisfactory description 
for both of these FF's. 

The investigation performed in this work has allowed us to describe the form factors 
F^yiQ^) (for g2 > 1.64 GeV^) and F^(Q2) (for > i.ig GeV^) in the context of the 
same theoretical scheme and by means of the same DA's. We have achieved a quite satisfac- 
tory agreement with the available experimental data. Theoretical computations have been 
carried out using the leading order [for F^{Q'^)] and the NLO [for F^^i^Q"^)] expressions for 
the hard-scattering amplitudes of the partonic subprocesses. An accuracy of our theoretical 
predictions may be improved by including into analyses the NLO and NNLO terms, respec- 
tively. These problems form directions for improving the developed theoretical framework 
and require separate detailed investigations. 
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